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We use 429 fb~ of e + e~ collision data collected at the T(4S) resonance with the BABAR detector 
to measure the radiative transition rate of b — > S7 with a sum of 38 exclusive final states. The 
inclusive branching fraction with a minimum photon energy of 1.9 GeV is found to be B(B — > 
X s j) = (3.29 ± 0.19 ± 0.48) x 10~ 4 where the first uncertainty is statistical and the second is 
systematic. We also measure the first and second moments of the photon energy spectrum and 
extract the best fit values for the heavy-quark parameters, mt and (1^, in the kinetic and shape 
function models. 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30. Er 

I. INTRODUCTION 

Flavor changing neutral current processes such as b — > 
S7, forbidden at the tree level in the standard model 
(SM), occur at leading order through radiative loop dia- 
grams. Since these diagrams are the dominant contribu- 
tions to this decay, the effects of many new physics (NP) 
scenarios, either enhancing or suppressing this transition 
rate by introducing new mediators within the loop, can 
be constrained by precision measurements of the total 
b — > S7 transition rate [IH5] . 

In the context of the SM, the first order radiative pen- 
guin diagram for the b — > 57 transition has a W boson 
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and f j or u quark in the loop. The SM calculation for 
the corresponding B meson branching fraction has been 
performed at next-to-next-to-leading order in the pertur- 
bative term, yielding B(B -> X s j) = (3.15±0.23) x 1(T 4 
for a photon energy of E 7 > 1.6 GeV, measured in the 
B meson rest frame u\. Experiments perform this 
measurement at higher minimum photon energies, gen- 
erally between 1.7 and 2.0 GeV, to limit the background 
from other B sources. The results are then extrapo- 
lated to the lower energy cutoff, _E 7 > 1.6 GeV, based 
on different photon spectrum shape functions. The cur- 
rent world average is in good agreement with the SM 
calculation, and is measured to be B(B — > Y s 7) = 
(3.55 ± 0.24 ± 0.09) x 10" 4 , for E 1 > 1.6 GeV 0. The 
second uncertainty is due to the photon spectrum shape 
function used to extrapolate to the 1.6 GeV photon en- 
ergy cutoff. 

The photon energy spectrum is also of interest, as it 
gives insight into the momentum distribution function 
of the quark inside the B meson. Precise knowledge of 
the function is useful in determining \V u b\ from inclusive 
semileptonic B — > X u lv measurements [9Tll3j. We fit the 
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measured spectrum to two classes of models: the "shape 
function" scheme [13] and the "kinetic" scheme [H] . The 
photon energy spectra predicted by these models are pa- 
rameterized to find the best values for the heavy quark 
effective theory (HQET) parameters, m& and /Lt^ [ID]. 

Our measurement uses a "sum of exclusives" approach, 
in which we reconstruct the final state of the s quark 
hadronic system, X s , in 38 different modes. For this ar- 
ticle we update a former BABAR analysis [TS] with about 
five times the integrated luminosity of the previous mea- 
surement, as well as an improved analysis procedure. By 
reconstructing the X s system, we access the photon en- 
ergy through: 



where E^ is the energy of the transition photon in the 
B rest frame, mj is the mass of the B meson, and mx s 
is the invariant mass of the X s hadronic system. Mea- 
suring mx, , with a resolution of around 5 MeV/c 2 , gives 
better resolution on E 7 than measuring the transition 
photon directly. We are also able to measure the en- 
ergy of the transition photon in the rest frame of the 
B meson rather than correcting for the boost of the 
B meson with respect to the center of mass (CM) as 
is required for a direct measurement of the transition 
photon. We perform this measurement over the range 
0.6 < mx 3 < 2.8GeV/c 2 in 14 bins with a width of 
100MeV/c 2S for m Xs < 2.0GeV/c 2 , and 4 bins with a 
width of 200MeV/c 2 for m Xs > 2.0GeV/c 2 . To eval- 
uate a total branching fraction for B(B — > X s j) with 
E 1 > 1.9 GeV, we sum the partial branching fractions 
from each mx 3 bin. This minimizes our dependence on 
the underlying photon spectrum structure, and is a de- 
parture from our previous procedure |15j . which com- 
bined the entire range 0.6 < mx s < 2.8GeV/c 2 and 
used a single fit to the signal yield to determine the total 
branching fraction. 



II. DETECTOR AND DATA 

Our results are based on the entire T(AS) dataset 
collected with the BABAR detector 16J at the PEP-II 
asymmetric-energy B factory at the SLAC National Ac- 
celerator Laboratory. The data sample has an inte- 
grated luminosity of 429 fb _1 collected at the T(4S) res- 
onance, with a CM energy y/s — 10.58 GeV, and contains 
471 x 10 6 BB pairs. We refer to this sample as the "on- 
peak" sample. An "off-peak" sample with an integrated 
luminosity of 44.8 fb^ 1 was recorded about 40 MeV be- 
low the T(45) resonance, and is used for the study of 
backgrounds consisting of e + e~ production of light qq (q 
=u, , s, ) or t + t~ states. 

The BABAR detector is described in detail in [16]. 
Charged-particle momenta are measured by the combina- 
tion of a silicon vertex tracker (SVT), consisting of five 



layers of double-sided silicon strip detectors, and a 40- 
layer central drift chamber (DCH) having a combination 
of axial and stereo wires. 

Charged-particle identification is provided by the com- 
bination of the average energy loss (dE/dx) measured in 
the tracking devices and the Cherenkov-radiation infor- 
mation measured by an internally reflecting ring-imaging 
Cherenkov detector (DIRC). 

Photon and electron energies are measured by a 
CsI(Tl) electromagnetic calorimeter (EMC). The SVT, 
DCH, DIRC, and EMC operate inside of a 1.5 T magnet. 
Charged n/fi separation is done using the instrumented 
flux return of the magnetic field, originally instrumented 
with resistive plate chambers [T5] and later with limited 
streamer tubes [T7]. 



III. SIGNAL AND BACKGROUND 
SIMULATION 

To avoid experimental biases, we use Monte Carlo 
(MC) simulations to model both the expected signal and 
background events and to define selection criteria before 
looking at the data. We have produced MC samples for 
e + e~ — > qq (q =u,,s,) and e + e~ — > t + t~ events, each 

at two times the on-peak luminosity, as well as BB MC 
events, excluding decays of the B meson to an X s "/ final 
state, at three times the on-peak luminosity. We also 
consider "cross-feed" backgrounds. We define cross-feed 
as signal events in which we wrongly reconstruct the B 
candidate. This occurs either because the X s final state 
is not one of the 38 reconstructed modes, not all of the 
particles in the true final state are detected, or the pro- 
cedure for selecting the correctly reconstructed B from 
several potential B candidates fails in some cases. 

Two types of signal MC events are generated, one for 
the if* (892) region (mx 3 < l.lGeV/c 2 ) in which the 
b — > ,S7 transition proceeds exclusively through B — > 
if* (892)7, and one for the region above the if* (892) 
resonance (1.1 < rnx s < 2.8 GeV/ c 2 , the upper bound 
being the limit of our ability to adequately reject B back- 
grounds). While there are several known X s resonances 
above the K* (892), and evidence for even more [IS], these 
resonances are broad and overlapping. We therefore take 
only the K* (892) resonance explicitly into account when 
simulating the signal events, as recommended by [19 . 

The quarks in inclusive region signal MC events shower 
using the "phase-space hadronization model," as opposed 
to the well known "Lund string model," with our de- 
fault JETSET [20] settings. The most important JET- 
SET parameters that influence the fragmentation of the 
X s system in this inclusive region are the probabilities 
of forming a spin-1 state for the s quark or uj quarks 
(the corresponding JETSET parameters are parj(12) 
and parj (11)). These probabilities are set to 0.60 and 
0.40, respectively. 

We generate the inclusive signal MC events with a flat 
photon spectrum with bounds corresponding to the mx. 



7 



boundaries, which we then reweight to match whichever 
spectrum model we choose. We do not take any ex- 
plicit photon model into account when evaluating sig- 
nal efficiency within a given X s mass bin. However, to 
evaluate the optimal background-rejection requirements, 
we do need to specify the expected shape of the spec- 
trum. For this, we use the model settings for the ki- 
netic scheme models [2] found to be consistent with 
the previous BaBAR sum of exclusive analysis (mj, = 
4.65GeV/c 2 , n% = 0.20 GeV 2 ) [H]. 

GEANT4 J5TJ is used to model the response of the 
detector for all MC samples. Time-dependent detector 
inefficiencies, monitored during data taking, are also in- 
cluded. 



IV. B MESON RECONSTRUCTION AND 
BACKGROUND REJECTION 

We reconstruct the B meson in one of 38 final states 
of the X s plus a high energy photon, as listed in Ta- 
ble [I] [55]. These modes consist of one or three kaons, at 
most one ?y, and at most four pions, of which no more 
than two can be neutral pions. The method of parti- 
cle identification (PID) has improved over the run of the 
experiment. In particular for charged K identification, 
we use a multi-class classifier procedure of error correct- 
ing output code (ECOC) [53]. The kaon identification 
efficiency is roughly 90% for the momentum range con- 
sidered for this analysis. 

TABLE I. The 38 X s decay modes used for B meson recon- 
struction in this analysis. 

Mode Num. Final State Mode Num. Final State 



1 




20 


K^TT + TT 7T + 7T 


2 


21 


K + iv + ir~ iv~ it 


3 




22 




4 


K s n 


23 


K+ V 


5 




24 


Kir, 


6 




25 


K° s riK + 


7 




26 


K + V n° 


8 




27 


K + r]-K~ 


9 


K+TT-TY 


28 




10 


ISO 
K S TT TV 


29 


K + rjTT + TV 


11 




30 


K° sV 7V + 7V° 


12 


K + 7T + 7T~~ 7T° 


31 


KItjty+ty- 


13 


K S TT TT 7T 


32 


K + r)-K~TT° 


14 




33 


K+K'K + 


15 




34 


K + K~K° S 


16 


A-+7r-7r 7r° 


35 


K + K-K° s k + 


17 


K + n + TV~ 7T + TV~ 


36 


K+K-K+n 


18 




37 


K + K~ K + n~ 


19 




38 


K + K-K° s n° 



The K® mesons are reconstructed as K® — > tt + tt can- 
didates with an invariant tt + tt~ mass within 9 MeV/c 2 of 
the nominal mass [IS], a flight distance greater than 



0.2 cm from the primary event vertex, and a flight sig- 
nificance (measurement of flight distance divided by the 
uncertainty on the measured distance) greater than 3. 
We do not include K® mesons or — > 7r°7r° decays in 
our reconstructed final states. 

Charged K candidates are identified based on the 
ECOC algorithms [53], which use information from the 
tracking system, the DIRC, and the EMC to identify par- 
ticle species using multivariate classifiers. All remaining 
charged tracks are assumed to originate from charged pi- 
ons. 

The 7T° and r\ candidates are reconstructed from pho- 
ton candidates with an energy greater than 60MeV as 
measured in the laboratory frame, and must have an in- 
variant mass between 115 and 150MeV/c 2 for the 7r°, and 
470 and 620MeV/c 2 for the r\. We also require a mini- 
mum momentum p^o „ > 200 MeV/c in the lab frame. Al- 
though we do not explicitly reconstruct the r, —> 7r + 7r~7r° 
decay mode, this mode is implicitly included in the final 
states if there is at most one other pion in the event. 
We combine these charged and neutral particles to form 
different X s candidates in the event. 

We require that an event contain at least one photon 
candidate with 1.6 < E* < 3.0 GeV (where "*" hence- 
forth indicates variables measured in the CM), which is 
consistent with the signal photon of the decay b — > sj. 
The distance to the closest cluster in the EMC is required 
to be greater than 25 cm from this signal photon clus- 
ter. We also require the angle between the signal photon 
candidate and the thrust axis of the rest of the event 
to satisfy |cos6^ 7 | < 0.85, and the ratio of event shape 
angular moments to satisfy L12/L10 < 0.46 [53] (the sig- 
nal peaks at slightly lower values than the background). 
These two preliminary requirements on the event topol- 
ogy are especially effective at reducing the large amount 
of more jet-like light qq backgrounds, and together de- 
crease this background source by about 50% (while only 
removing 10% of signal). 

We combine the X s candidates and the signal pho- 
ton candidates to form B candidates in the event. 
We define the beam-energy substituted mass, mgs = 
a/( v / s/2) 2 - (pb) 2 , and require Toes > 5.24 GeV/ c 2 . We 
also require the difference between the expected B energy 
and the reconstructed B energy, \AE\ = \E* B — \/s/2\, 
to satisfy \ AE\ < 0.15 GeV. For these quantities, p* B and 
E* B are the momentum and energy of the reconstructed 
B meson in the CM system. 

With these loose preliminary requirements in place, 
each event still typically has several B meson candidates. 
We construct a random forest classifier [5^] (a signal se- 
lecting classifier, SSC) to find the best candidate in an 
event. This classifier is built using the variables AE/ge 
(where ue is the uncertainty on the total energy of the 
reconstructed B), the thrust of the reconstructed B, the 
7T° momentum in the CM (if the candidate has a tt°), the 
invariant mass of the X s candidate, and the zeroth and 
fifth Fox- Wolfram moments of the event (5^] . We choose 
to include the fifth Fox- Wolfram moment because our 



MC simulation indicates that this variable improves the 
performance of our classifier. The selected B candidate 
in an event is the candidate with the highest response to 
this classifier. We find that applying this classifier to se- 
lect the best candidate, after placing a loose requirement 
on | AS |, rather than selecting the candidate with the 
smallest lA^I, improves the signal efficiency by a factor 
of about two. We also find that placing a requirement 
on the SSC response is effective at further removing B 
backgrounds. 

To further reduce the background from events in which 
a photon from a high energy ir° decay is mistaken as 
the signal photon candidate, we construct a dedicated 
7r° veto using a random forest classifier [35] . If the signal 
photon candidate in an event can be combined with any 
other photon to form a candidate with an invariant mass 
in the range 115 < m 77 < 150MeV/c 2 , we evaluate the 
7T° veto classifier response based on the invariant mass 
of the two photons and the energy of the lower energy 
photon. The response of the ir° veto classifier is used as 
input to a more general background rejecting classifier 
(BRC). 

The BRC is constructed to remove continuum (lighter 
qq) backgrounds. To construct this classifier, we use in- 



for the combinatorial components: 



formation from the it veto, I cos 9, 



T 7 | 



COS 



between the thrust axis of the B and the thrust axis of 
the rest of the event), |cos(9^| (the CM polar angle of the 
B flight direction), the zeroth, first, and second angular 
moments |24| computed along the signal photon candi- 
date's axis as well as the ratio ^12/^10 (which exhibits 
slightly different signal and background shapes), and the 
10° momentum flow-cones around the B flight-direction. 

To effectively remove background while maintaining 
signal efficiency, we evaluate optimal requirements for the 
responses of the BRC and SSC in four mass regions, [0.6- 
1.1], [1.1-2.0], [2.0-2.4], and [2.4-2.8] GeV/c 2 , optimizing 
the figure of merit S/y/S + B, where 5* is the expected 
signal yield and B is the expected background yield eval- 
uated using MC simulation. 



V. SIGNAL YIELD EXTRACTION 

We extract the signal yield by performing fits to the 
Toes distribution in each bin of mx B ■ The signal distri- 
bution is described by a Crystal Ball function (CB) |27j : 



/(m E s) = eS 



' ncB\ ncB 



/(m E s) 




< a, 



> a, 



(2) 



where m and a are the peak position and width, respec- 
tively, and the parameters a and hqb take account of the 
non-Gaussian tail. This distribution takes into account 
the asymmetry of the Toes distribution for these events. 
The backgrounds are described by ARGUS functions [25] 



f(m E s) = m ES 1 - 



/ rn E s V 
Km) 



x e v 



(3) 



where m is the end point, and c is the slope, and Novosi- 
birsk functions [29 for both the peaking BB contribution 
and peaking cross- feed contribution ( "peaking" meaning 
apparently resonant behavior similar to the signal distri- 
bution in toes)- 

The signal CB distribution is parameterized based on 
a fit to correctly reconstructed signal MC events over 
the full hadronic mass range, 0.6-2.8 GeV/c 2 , as we find 
little X s mass dependence of the signal shape parameters. 
The CB parameters take the values a=1.12, too=5.28 
GeV/c 2 , <t=2.84 MeV/c 2 , and n CB = 145 for every mass 
bin. In Sec. lVIII we evaluate the uncertainties indroduced 
by fixing the CB shape parameters. 

The cross-feed shape has both a peaking component 
and a combinatoric tail. The peaking component is de- 
scribed by a Novosibirsk function, parameterized over 
five different mass regions, [0.6-1.1], [1.1-1.5], [1.5-2.0], 
0£| (the angle [2.0-2.4], and [2.4-2.8] GeV/c 2 , based on MC distribu- 



tions over these regions. The combinatoric cross-feed tail 
is described by an ARGUS function with the slope c fit 
to the MC events in each mass bin, and fixed in the fits 
to data. We fix the fraction of peaking cross-feed MC 
events, the fraction of signal to signal+cross-feed events, 
and the shapes of the cross-feed Novosibirsk and ARGUS 
functions, in each bin of mx s , based on the MC events. 
We allow the total signal+cross-feed yield to float in each 
mass bin in the fits to data. 

A second ARGUS function is used to parameterize the 
combinatoric background from continuum and other BB 
sources. We fix the end point to of the ARGUS function 
to the kinematic limit (5.29GeV/c 2 ) of the toes variable 
and allow the yield to float. 

The BB background also has a peaking component, 
which becomes more significant at higher X s mass, is 
also described by a Novosibirsk function, and is param- 
eterized over three mass ranges, [0.6-2.0], [2.0-2.4], and 
[2.4-2.8] GeV/c 2 . We fix the total number and shape of 
the peaking BB events based on a fit to the BB MC 
events in each mass bin. 

We perform a minimum y 2 fit to the toes distribution 
in each bin of mx 3 , allowing the slope of the combinatoric 
ARGUS and the fractional yield of signal+cross-feed to 
float (the complementary fractional yield, once the peak- 
ing BB is accounted for, reflects the normalization of 
the combinatoric ARGUS function) . Figure [l] shows an 
example for mx s bin 1.4-1.5 GeV/c 2 . We fix all other 
shape parameters, and evaluate systematic uncertainties 



associated with fixing these parameters in Sec. VII We 
perform MC simulations ("toy studies") to ensure that 
we do not introduce any biases due to the fitting proce- 
dure. 
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FIG. 1. The fit in mass bin 1.4 < m Xs < 1.5GeV/c 2 to 
(a) signal MC events, (b) cross-feed MC events, (c) peaking 
BB MC events, and (d) the data. The signal (thick dashed), 
cross-feed (two dot-dashed curves, one ARGUS function and 
one Novosibirsk function), peaking BB (dotted), and com- 
binatoric background (thin dashed) component functions are 
shown. 



VI. X s FRAGMENTATION AND MISSING 
FRACTION 

The fragmentation of the hadronic system in the in- 
clusive region, 1.1 < mx s < 2.8GeV/c 2 , is modeled with 
JETSET with a phase-space hadronization model. The 
differences between fragmentation in the MC sample and 
in the data influence the measurement in two ways. First, 
since the efficiencies for the 38 modes are not the same, 
an incorrect modeling of their relative fractions will lead 
to an incorrect expected total efficiency for reconstruct- 
ing the 38 final states (£38 )■ Second, the simulation of the 
fragmentation process can introduce incorrect estimates 
of the fraction of the total inclusive b — > sj transition 
rate reflected by the 38 modes (e^d). The fraction of 
final states in each of the mass bins that is not included 
in our 38 modes is referred to as the "missing fraction" , 
and is equivalent to 1 — £i nc i- 

We are able to evaluate and correct 63s for the first ef- 
fect, and use these results to estimate the uncertainty on 
the second effect, our uncertainty on ei ac \, by performing 
a fragmentation study comparing the frequency of groups 
of modes in the MC sample to the data. For this study, 
we compare the frequency of ten groups of modes, each 
containing two to ten final states, in the MC sample to 
the frequency for these groups found in the data. We per- 
form this study in four different mass regions, [1.1-1.5], 
[1.5-2.0], [2.0-2.4], and [2.4-2.8] GeV/c 2 . 

The procedure for the study involves reweighting the 
relative contribution of each of the groups of modes in 
our MC based on the relative amount found in the data. 
The efficacy of the procedure is checked on MC events 
by ensuring we can find the £38 in each mass bin for the 
Lund string model when starting with the default phase- 
space hadronization model [50], as well as find the £33 in 
each mass bin for the phase-space hadronization model 
when starting with the Lund string model. The different 
groups of modes we use to compare data and the MC 
samples, along with the results of the comparisons in 
each mass bin, are given in Table [TTJ, obtained with the 
default phase-space hadronization model as the starting 
point. 

To perform this study, we combine the mass bins into 
the four mass regions, and fit the signal+crossfeed con- 
tribution for each subset of modes in each mass region in 
the data. We then use the ratio of the yield of each subset 
found in data to the amount found in the MC sample to 
reweight the MC sample to better reflect the data in the 
mass region. Wc use the statistical uncertainty in fitting 
each subset in data as the uncertainty on the ratio. 

After correcting the signal and cross-feed MC events 
based on these comparisons, we evaluate the value of £38 
for each mass bin, reported in Table \LU\ For the inclusive 
region, the uncertainty on £38 is calculated using the un- 
certainties in the fragmentation corrections, as described 
later in Sec. VII Since the fragmentation in the K* (892) 
region is considered well modeled, we do not perform a 
fragmentation correction on these mass bins. 
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TABLE II. The subsets of modes and the ratio of yields found in each mx 3 region when comparing the data to the MC events. 
The error is statistical only. 



Data 




Definition 




Modes used 


subset 










1 




2 bodies without 7r u 




1,3 


2 




2 bodies with 1 it 




2,4 


3 




3 bodies without ty° 




5,8 


4 




3 bodies with 1 it 




6,9 


r. 





4 bodies without 7r^ 




11 1/1 


6 




4 bodies witn 1 n 




12,15 


n 
1 

8 




o/4 bodies with 2 7r 
5 bodies with 0-2 7r 




7,10,lo,16 
17-22 


n 




J] — y ryry 




16-61 


i n 
1U 




olv modes 




oo oo 
00-08 


Data 
subset 


1.1 < mx B < l.5GeV/(? 
(data/MC) 


1.5 < m Xs < 2.0GeV/c 2 
(data/MC) 


2.0 < m Xs < 2.4GeV/c 2 
(data/MC) 


2.4 < m Xs < 2.8 GeV/c 2 
(data/MC) 


1 


0.65 ± 0.03 


0.38 ± 0.03 


0.05±0.05 


0.18± 0.13 


2 


0.53 ± 0.05 


0.28 ± 0.06 


0.32 ± 0.12 


U.IO _ .15 


3 


1 9(1 + OS 

_L.ZiU _L_ U.Ud 


1 01 + 04 


79 + 1 1 


95+n 25 


4 


1.70 ± 0.05 


1.03 ± 0.06 


0.33 ± 0.13 




5 


0.34 ± 0.08 


1.34 ± 0.10 


1.12 ± 0.23 


2.29 ± 0.74 


6 


1.24 ± 0.13 


1.16 ± 0.11 


1.28 ± 0.27 


in+ - 39 

u - lu -0.10 


7 


0.56 ± 0.19 


1.37 ± 0.30 


0.83 ± 0.53 


2.06 ± 1.64 


8 


1 on+ 105 

x - uu -1.00 


0.57 ± 0.16 


0.74 ± 0.28 


o.29iJ;S 


9 


0.94 ± 0.15 


1.70 ± 0.20 


2.47 ± 0.50 


1 09+ 103 


10 


0.00 ± 0.00 


0.62 ± 0.11 


0.74 ± 0.31 


0-83_ ' 83 



TABLE III. The value of 63s before and after the fragmenta- 
tion corrections are performed on the inclusive region. The 
uncertainty on the corrected value in the inclusive region re- 
flects the uncertainty of the fits to the data. 





e 3 8 original 


£38 final 


(GeV/c 2 ) 


(%) 


(%) 


0.6-0.7 


15.0 


15.0 


0.7-0.8 


16.5 


16.5 


0.8-0.9 


17.3 


17.3 


0.9-1.0 


18.3 


18.3 


1.0-1.1 


16.0 


16.0 


1.1-1.2 


11.5 


10.4±0.4 


1.2-1.3 


11.6 


10.6±0.3 


1.3-1.4 


10.7 


9.9±0.3 


1.4-1.5 


9.5 


8.9±0.5 


1.5-1.6 


8.4 


7.5±0.5 


1.6-1.7 


7.2 


6.5±0.4 


1.7-1.8 


5.5 


5.0±0.4 


1.8-1.9 


4.5 


4.2±0.4 


1.9-2.0 


3.3 


3.0±0.4 


2.0-2.2 


4.0 


3.2±0.4 


2.2-2.4 


3.1 


2.4±0.4 


2.4-2.6 


2.3 


1.9±0.7 


2.6-2.8 


2.3 


2.1±0.9 



We base the uncertainty on the fraction of the inclu- 
sive b — > S7 transitions measured by the 38 final states, 
£i nc i, on the range of values predicted by competing frag- 
mentation models in the MC samples. We consider many 



settings of JETSET using both the default phase-space 
and the Lund string hadronization mechanism, as well 
as a thermodynamical model 30j . Other models in JET- 
SET (Field-Feynman model j31] of the showering quark 
system, etc.) are found to yield results consistent with 
the Lund string model, and are not further considered. 

As mentioned above, we identify the probabilities for 
forming a spin-1 hadron with the s quark or it/quarks to 
be the JETSET parameters that have the largest impact 
on the breakdown of final states. We try many settings 
for these parameters in both the phase-space hadroniza- 
tion mechanism and the Lund string model mechanism in 
JETSET. By varying the spin-1 probabilities and using 
both of these fragmentation mechanisms, we are able to 
identify a range of models that, taken together, account 
for the breakdown of final states found in the data in the 
fragmentation study (Table [TT| . We vary the probability 
for forming a spin-1 hadron with the s quark between 
zero and one, and with the it/quark between 0.2 and 0.8. 
When comparing to our default MC settings, the mod- 
els we consider predict both higher and lower ratios than 
those found in the data, but no single model matches 
every ratio in every mass region. 

We also find that no single mechanism or JETSET set- 
ting perfectly reproduces the fragmentation in the data; 
however the models chosen bound the data. The fact that 
spin-1 probability settings need to be varied to account 
for data and MC differences is expected, as a variety 
of resonances exist in the inclusive region. The maxi- 
mum, minimum, and default values for ei nc i that we find 
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are reported in Table IV We account for what is seen 



in data in the fragmentation study through a variety of 
settings of both the Lund string mechanism and phase- 
space hadronization mechanism, and therefore base our 
uncertainty on ei nc i on these same models. The statistics 
model and the default JETSET settings predict values for 
ej nc i that lie between those predicted by other settings of 
JETSET that we try. As we find that no model exactly 
describes the fragmentation we observe in the data, but 
together the models considered bound the data, we count 
each model as equally probable, and take the systematic 
uncertainty on the correct value for £i nc i as the differ- 
ence between the maximum and minimum values of £i nc i 
relative to the default MC value, and divide by V12, re- 
flecting the standard deviation of a uniform distribution. 



TABLE IV. The minimum, maximum, and default values of 
ei nc i found for the range of models that account for the dif- 
ferences seen between the default MC events and data in the 
inclusive region. We include the 7f*(892) region default val- 
ues as well, though these mass bins are not modeled by the 
inclusive MC sample. 



mx s Minimum Maximum Default 

(GeV/c 2 ) e inc i einci e in ci 

0.6-0.7 0.75 

0.7-0.8 0.74 

0.8-0.9 0.74 

0.9-1.0 0.75 

1.0- 1.1 0.74 

1.1- 1.2 0.71 0.74 0.73 

1.2- 1.3 0.71 0.74 0.72 

1.3- 1.4 0.70 0.74 0.72 

1.4- 1.5 0.69 0.73 0.71 

1.5- 1.6 0.66 0.73 0.68 

1.6- 1.7 0.59 0.72 0.66 

1.7- 1.8 0.57 0.72 0.63 

1.8- 1.9 0.52 0.71 0.59 

1.9- 2.0 0.47 0.68 0.54 
2.0-2.2 0.41 0.64 0.48 
2.2-2.4 0.33 0.60 0.39 
2.4-2.6 0.27 0.56 0.31 
2.6-2.8 0.23 0.51 0.25 



VII. SYSTEMATIC UNCERTAINTIES 

We present the X s mass-bin-dependent uncertainties 
in Table |Vj The uncertainty on the total number of B 
mesons produced at BaBAR is evaluated at 1.1%. 

The uncertainty on the efficiency of the requirements 
on the two multivariate classifiers are evaluated in signal- 
like data sidebands, regions in parameter space similar 
to, but not overlapping with, the signal region, by com- 
paring the efficiency of the requirements on MC events 
and the efficiency of these requirements on data. We 
define our sidebands as the inverse of the requirements 
we place on the classifiers. Therefore if we require the 



SSC response to be greater than 0.5 for a mass region, 
we evaluate the BRC uncertainty in the SSC sideband 
defined by requiring an SSC response less than 0.5 (and 
similarly for evaluating the SSC uncertainty in the BRC- 
defined sideband). The relative difference between the 
two efficiencies is taken as the systematic uncertainty. 
The sideband produced by taking the inverse of the re- 
quirements on the SSC is used to evaluate the uncertainty 
on the requirements on the BRC. 

To evaluate the uncertainty on the SSC requirement, 
the events that are identified by the 7r°-veto classifier to 
contain a 7r° candidate are used with the further require- 
ment rriES > 5.27GeV/c 2 . This gives a more signal-like 
sample of events that have a high energy 7r° in place of 
the signal transition photon. The efficiency of the SSC re- 
quirement is compared between data and the MC events 
with the use of this sideband. 

To evaluate fitting uncertainties related to fixing many 
of the parameters in the signal and cross-feed PDFs, we 
use the iV* (892) region (rax, < LI GeV/c 2 ) to determine 
reasonable shifts in these parameters. We assign the sys- 
tematic uncertainty as the change in signal yield in the fit 
to data when we use the shifted shape parameters. For 
the uncertainty on the fraction of signal to signal+cross- 
feed, which is also fixed in the fit to data, we fix the total 
yield and slope of this ARGUS function (these are the 
two parameters that we float in the fits to data) and al- 
low this fraction to float in each mass bin. We take the 
change in signal yield when we fix the signal fraction to 
this new value as the systematic uncertainty. 

To evaluate the uncertainty on the peaking BB back- 
ground PDF shape, we use the change in signal yield 
when changing the parameter values by the uncertainty 
in the fits to MC events. 

The uncertainty on the number of peaking BB events, 
generally the largest source of BB fitting error in Ta- 
bic |vj is again evaluated based on the 7r°-veto sideband. 
In this sideband, we evaluate the BB MC predictions for 
the number of peaking events and compare this to the 
number of peaking BB events we find in data. We find 
these values to agree within one standard deviation for 
the three mass regions over which we have parameter- 
ized the peaking BB Novosibirsk function (see Sec. |v|). 
We determine the mass-region-dependent uncertainty on 
the measurement of peaking BB yield in the 7r°-veto side- 
band in data. We use this uncertainty added in quadra- 
ture with the uncertainty from the fits to the BB MC 
sample as the uncertainty on the number of peaking BB 
events in each mass bin. Unlike the other systematic 
uncertainties, which are multiplicative in nature, this un- 
certainty is additive since we are subtracting out peaking 
BB events we would otherwise fit as signal+cross-feed in 
the fits to data. 

The detector response uncertainties associated with 
PID, photon detection both from the transition photon 
and from ir°/r] decay, and tracking of charged particles 
are approximately 2.5-2.9% in every mass bin. 

The uncertainty on 63$ from the fragmentation study is 
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taken from the change in 63s when modifying the weights 
given in Table [n] by the uncertainty on these values in- 
dividually. We also account for the differences in statis- 
tics between the mass regions over which these uncer- 
tainties were determined and the individual mass bins. 
Since our fragmentation study procedure groups bins to- 
gether before evaluating appropriate weights, the weights 
we identify tend to reflect the bins with higher numbers of 
events, and the uncertainty on the bins with fewer events 
needs to be increased. We therefore increase the uncer- 
tainty in each rrix s bin by a factor of J iV rcg i on / y 7 A^m , 
where iV rog i on (A^in) refers to the number of events in 
the region (bin). This correction ensures that if an mx 3 
bin has few events compared to its corresponding re- 
gion, then the uncertainty for this bin will be larger. 
The total fragmentation uncertainty is found by sum- 
ming in quadrature the changes for each of the ten subset 
amounts. Where asymmetric uncertainties are reported 
in Table |IT} we take the average change in 63% when fluc- 
tuating the weights by the indicated amounts. For the 
mass bin 1.0 < mx B < 1-lGeV/c 2 , it is unknown if the 
fragmentation in the data is modeled more effectively by 
the AT* (892) MC sample or the inclusive MC sample. We 
take the average of the two predictions to be the value 
for €38, and the uncertainty is the difference divided by 
VT2, consistent with the standard deviation of a uniform 
distribution. 



The uncertainty on the missing fraction was covered in 
Sec. [VI] for the inclusive region. The competing fragmen- 
tation models give an uncertainty on the missing fraction 
from 1.3 to 32.7%, getting larger at higher mass. For the 
if* (892) region, we take the uncertainty to be the differ- 
ence between the default if* (892) MC prediction for the 
missing fraction, and the hypothesis of exclusively miss- 
ing K® final states, which would be a missing fraction of 
25% for this region. 



We take each of these systematic uncertainties to be 
uncorrelated within an mx s bin. However, there are 
correlations in the errors between the mass bins. The 
BB counting, classifier requirements, non-BB fitting for 
signal and cross-feed PDF shape, and detector response 
systematic uncertainties are taken to be completely cor- 
related between all mass bins. As we parameterize the 
peaking BB Novosibirsk function in three different re- 
gions, we evaluate the uncertainties over the same re- 
gions, taking the uncertainties to be independent from 
one region to the other (indicated by the horizontal lines 
in Table [\7|. Similarly, the fragmentation uncertainty and 
missing fraction uncertainty are evaluated using different 
samples and strategies in different mass regions; we take 
the uncertainty on these mass regions to be uncorrelated 
with one another, but completely correlated between the 
mass bins within a mass region. 



VIII. BRANCHING FRACTIONS 

We measure the signal yield in 100 MeV/c 2 wide bins 
of the X s mass over the range 0.6 < mx„ < 2.0GeV/c 2 , 
and 200 MeV/c 2 wide bins over the mass range 2.0 < 
mx s < 2.8 GeV/c 2 . We report the measured signal yield 
in Table VI where we have included the \ 2 per degree 
of freedom (dof) from the fits. 

We use the efficiencies reported in Tables |Hl| and |IV| to 
derive the total number of b — > sj events, Nb^ S7 , based 
on the yields, A/yi e id, reported in Table VI according to: 



N, 



b— >S7 



Afield 
£38£incl 



(4) 



The partial branching fraction (PBF) for each mass bin 
is reported in Table |VII| In this table, we also report the 
total branching fraction, with a minimum photon energy 
of E-f > 1.9 GeV, reflecting the sum of the 18 bins: 



B(B -> X sl ) = (3.29 ± 0.19 ± 0.48) x 10" 



(5) 



where the first uncertainty is statistical and the second 
is systematic. This result is consistent with the previ- 
ous BaBAR sum of exclusives results of B(B — > X s j) — 
(3.27 ± 0.18±j^±^) x 10" 4 Q2], where the first un- 
certainty is statistical, the second systematic, and the 
third from theory. The total statistical uncertainty on 
our result reflects the sum in quadrature of the statisti- 
cal uncertainty of the 18 uncorrelated statistical uncer- 
tainties in the mass bins. This method ensures reduced 
spectrum model dependence when quoting a branching 
fraction. An alternate method of measuring the transi- 
tion rate based on larger mass bins yields similar results. 
This alternative method is similar to the method used in 
the previous analysis |15j . in which one measurement of 
the signal yield over the entire mass range was used to de- 
termine the total transition rate. However, that method 
introduces additional model dependence due to the un- 
certainty in the spectrum shape and we instead decide 
to take the total transition rate as the sum of the transi- 
tion rates in each of the mx s bins. The total systematic 
uncertainty reported in our study takes the correlations, 
indicated in Table |V} into account. The correlation coef- 
ficients between the total uncertainties in each bin are in- 
cluded in the Appendix. The partial branching fractions 
per 100 MeV/c 2 in X s mass are illustrated in Fig. [2J with 
the previous BABAR sum-of-exclusive results also shown. 



IX. FITS TO SPECTRUM MODELS AND 
MOMENTS 



Since we measure the hadronic mass spectrum in bins 
of 100 or 200 MeV/c 2 , we are able to fit directly different 
models of this spectrum to obtain the best-fit values of 
different HQET parameters. We choose to fit two such 
classes of models: the kinetic model using an exponential 
distribution function [14] , and the shape function model, 
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TABLE V. List of systematic uncertainties described in the text. These sub-component systematic uncertainties are assumed 
to be uncorrelated within a given mass bin and the total uncertainty reflects their addition in quadrature. All uncertainties 
are given in %. Many of these uncertainties are taken to be completely correlated over mx s regions, and we have indicated the 
correlated uncertainties with horizontal lines defining the regions. 

Mass bin BB Classifier Non-SB BB Detector Frag. Missing Total 



GeV/c 2 ) 


counting 


selection 


fitting 


fitting 


response 




fraction 




0.6-0.7 


1.1 


1.0 


14.9 


21.3 


2.5 




0.6 


26.2 


0.7-0.8 


1.1 


1.0 


2.7 


3.1 


2.6 




0.9 


5.1 


0.8-0.9 


1.1 


1.0 


1.7 


0.6 


2.6 




1.3 


3.8 


0.9-1.0 


1.1 


1.0 


1.7 


0.7 


2.7 




0.0 


3.6 


1.0-1.1 


1.1 


1.0 


5.1 


2.5 


2.7 


13.1 


0.9 


14.6 


1.1-1.2 


1.1 


0.7 


5.7 


0.9 


2.7 


3.9 


1.3 


7.7 


1.2-1.3 


1.1 


0.7 


4.7 


0.4 


2.7 


3.0 


1.3 


6.4 


1.3-1.4 


1.1 


0.7 


4.6 


0.3 


2.7 


3.0 


1.6 


6.4 


1.4-1.5 


1.1 


0.7 


4.7 


0.6 


2.7 


5.7 


1.8 


8.2 


1.5-1.6 


1.1 


0.7 


3.7 


1.5 


2.7 


6.1 


3.1 


8.5 


1.6-1.7 


1.1 


0.7 


4.3 


1.3 


2.7 


6.3 


5.9 


10.2 


1.7-1.8 


1.1 


0.7 


4.9 


1.5 


2.7 


7.9 


6.9 


12.1 


1.8-1.9 


1.1 


0.7 


3.4 


13.1 


2.7 


10.0 


9.6 


19.6 


1.9-2.0 


1.1 


0.7 


5.3 


4.2 


2.7 


13.4 


11.1 


18.9 


2.0-2.2 


1.1 


1.9 


4.5 


6.6 


2.9 


11.0 


13.9 


19.8 


2.2-2.4 


1.1 


1.9 


4.9 


22.0 


2.9 


18.4 


19.7 


35.3 


2.4-2.6 


1.1 


2.8 


4.7 


23.8 


2.8 


36.7 


26.8 


51.7 


2.6-2.8 


1.1 


2.8 


49.3 


154.1 


2.8 


45.7 


32.7 


171.3 



TABLE VI. Signal yields from fits to the on-peak data and 
corresponding x 2 /dof from the fits (the uncertainties are sta- 
tistical only). 



rax s 


-^Vyield 


Data fit 


(GeV/c 2 ) 


(events) 


X 2 /dof 


0.6-0.7 


5.9 


± 


12.2 


0.8 


0.7-0.8 


114.7 


± 


24.0 


0.9 


0.8-0.9 


2627.4 


± 


50.2 


1.0 


0.9-1.0 


2249.5 


± 


53.1 


0.9 


1.0-1.1 


380.4 


± 


36.1 


0.9 


1.1-1.2 


393.7 


± 


37.1 


0.8 


1.2-1.3 


1330.5 


± 


47.1 


0.6 


1.3-1.4 


1501.0 


± 


54.7 


1.0 


1.4-1.5 


1479.6 


± 


58.3 


1.0 


1.5-1.6 


1039.6 


± 


55.7 


0.9 


1.6-1.7 


929.1 


± 


56.7 


0.9 


1.7-1.8 


736.5 


± 


48.6 


1.2 


1.8-1.9 


585.8 


± 


50.8 


1.0 


1.9-2.0 


272.0 


± 


37.4 


1.0 


2.0-2.2 


684.4 


± 


68.2 


1.1 


2.2-2.4 


277.5 


± 


64.6 


1.0 


2.4-2.6 


159.7 


± 


54.4 


0.8 


2.6-2.8 


-34.4 


± 


62.0 


1.1 



also using an exponential distribution function [13]. The 
choice of distribution function is not expected to have a 
large impact on the values determined for the underly- 
ing HQET parameters for each model, but the parame- 
ters themselves arc not immediately comparable between 
models (for example, the models are evaluated at differ- 
ent energy scales). 



TABLE VII. The partial branching fractions in each mass bin 
reflecting branching fractions per 100 or 200 MeV/ c 2 , and the 
total branching fraction for b — > S7 with > 1.9 GeV. The 
uncertainties quoted are statistical and systematic. 



nil, Branching fraction 
(GeV/c 2 ) per 100 or 200 MeV/c 2 
(xKT 6 ) 



0.6-0.7 


0.1 


± 


0.1 


± 


0.0 


0.7-0.8 


1.0 


± 


0.2 


± 


0.1 


0.8-0.9 


21.8 


± 


0.4 


± 


0.8 


0.9-1.0 


17.4 


± 


0.4 


± 


0.6 


1.0-1.1 


3.4 


± 


0.3 


± 


0.5 


1.1-1.2 


5.5 


± 


0.5 


± 


0.4 


1.2-1.3 


18.4 


± 


0.7 


± 


1.2 


1.3-1.4 


22.5 


± 


0.8 


± 


1.5 


1.4-1.5 


24.9 


± 


1.0 


± 


2.0 


1.5-1.6 


21.5 


± 


1.2 


± 


1.8 


1.6-1.7 


23.0 


± 


1.4 


± 


2.3 


1.7-1.8 


24.6 


± 


1.6 


± 


3.0 


1.8-1.9 


25.4 


± 


2.2 


± 


5.0 


1.9-2.0 


17.9 


± 


2.5 


± 


3.4 


2.0-2.2 


24.0 


± 


2.4 


± 


4.7 


2.2-2.4 


16.2 


± 


3.8 


± 


5.7 


2.4-2.6 


14.1 


± 


4.8 


± 


7.3 


2.6-2.8 


-3.5 


± 


6.4 


± 


6.1 


0.6-2.8 


329 


± 


19 


± 


48 



In order to fit the measured spectrum to these mod- 
els, we need to take special account of the AT* (892) res- 
onance, as the models assume quark-hadron duality in 
their spectra. Consequently, the models smooth over this 
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FIG. 2. The partial branching fractions binned in (a) X s mass 
and (b) the corresponding bins, with the statistical and 
systematic uncertainties added in quadrature. The current 
results (solid lines) and former BABAR results [T5] (dashed 
lines) are shown. 



resonance. We fit a relativistic Breit-Wigner [32] (RBW) 
to the K* (892) MC sample at the generator level to ex- 
tract the parameters of this curve. Fits to the transi- 
tion point between the RBW curve of the if* (892) reso- 
nance and the remaining spectrum indicate a value close 
to mx s =l-17 GeV/c 2 , which we take to be the location 
of this transition. Furthermore, we require that the in- 
tegral of the RBW used to parameterize the if* (892) 
region (rax, < 1.17GeV/c 2 ) be equivalent to the in- 
tegral of this region in the spectrum models. For the 



hadronic mass bin containing the transition from the 
if* (892) resonance to the nonresonant-spectrum mod- 
els (1.1 < m Xs < 1.2 GeV/ c 2 ), we assign the value 
of the integral of the RBW up to the transition point 
(1.10 < rnx s < 1.17 GeV/ c 2 ) plus the integral of the spec- 
trum model from the transition point to the bin boundary 
(1.17 < m Xs < 1.20GeV/c 2 ). 

We perform a fit to the different spectrum models by 
minimizing the quantity 



E 



(PBF th - PBF exp ) 4 C'-/ (PBF th - PBF exp ) 7 



(6) 

where PBF t h and PBF exp are the PBF predicted by the 
spectrum model in the mass bin and the PBF we mea- 
sured in the mass bin, respectively. The matrix 
is the inverse of the matrix of correlation coefficients 
between the uncertainties on bins i and j, reported in 
Appendix [XJ having taken the correlated systematic un- 
certainties and uncorrelated statistical uncertainties into 
account. The <7j and <jj are the total uncertainties (statis- 
tical and systematic added in quadrature) on the branch- 
ing fractions determined for bins i and j. 

We find the best HQET parameter values based on 
the measured hadronic mass spectrum for two quantities 
for each model we fit. For the kinetic model, we fix the 
chromomagnetic operator (p%) to 0.35 GeV 2 , and the 
expectation values of Darwin (p% = 0.2 GeV 3 ) and spin- 
orbit (p 3 LS = -0.09 GeV 3 ) terms; we allow m b and /x 2 to 
take values between 4.45 and 4.75GeV/c 2 and 0.2 and 
0.7 GeV 2 , respectively. We have points on the m b -p1 
plane at which the spectrum has been evaluated exactly. 
These points are spaced every 0.05 GeV/ c 2 for m b and 
every 0.05 GeV 2 for /i 2 . We interpolate the spectrum 
mass bin predictions between these points using 



F(m b , pl) = A + Bx (m b - 4.45) + C x (pi - 0.2) + D x (m b - 4.45)(/4 - 0.2), 



(7) 



where we solve this equation for [A, B,C, D]. The val- 
ues 4.45 and 0.2 in Eq.Q are changed to the different 
values for which we have exact spectra provided. This 
strategy ensures continuity in the value of the spectrum 
predictions for each hadronic mass bin across the m b -p^ 
plane. 

The shape function models use two variables to param- 
eterize the spectrum, b and A [13] , that may be converted 
to values of m b and /i 2 , evaluated at a single energy scale 
of 1.5 GeV. Similar to the kinetic model fits, we interpo- 
late between points on the b-A plane at which we have 
exact spectrum predictions (for 2.0 < b < 5.0 in incre- 
ments of 0.25, and 0.4 < A < 0.9 GeV in increments of 
0.05 GeV). 



The best fit values for the HQET parameters are re- 
ported in Table |VIII| The uncertainty reflects the values 
at which the value of x 2 changes by one unit. The cor- 
responding best fit spectrum model and la error ellipses 
are shown in Fig. [3] (kinetic model) and|l] (shape function 
model) . 

TABLE VIII. The best fit HQET parameter values based on 
the measured mx s spectrum. 



Kinetic model |14| 


Shape function model 


m 


m h 4.568±H!e GeV/c 2 
l4 0.450 ±0.054 GeV 2 


4.579±r ^GeV/ C 2 
0.257l° ;^GeV 2 
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FIG. 3. The (a) One-o" region for the kinetic model param- 
eters based on the measured spectrum and the (b) best fit 
kinetic model compared to the measured PBFs. The error 
bars in (b) include the statistical and systematic errors added 
in quadrature. 



FIG. 4. The (a) One-<r region for the shape function model 
parameters based on the measured spectrum and the (b) best 
fit shape function model compared to the measured PBFs. 
The error bars in (b) include the statistical and systematic 
errors added in quadrature. 



We use the PBFs measured in each mass bin to cal- 
culate the mean and variance of the photon energy spec- 
trum. These quantities are spectrum-model independent, 
and may be used to constrain the parameters in other 
models. We evaluate the mean and variance for five dif- 
ferent minimum photon energies and report the values in 
Table llXl 



uncertainties on the mean and variance calculated at dif- 
ferent photon energy cutoffs. We report these values in 
Table [X] When determining the uncertainty on the means 
and variances, and evaluating the correlations between 
these uncertainties, we take into account the correlated 
systematic errors reported in Table |V| 



TABLE IX. The mean and variance of the photon energy spec- 
trum, calculated for five photon energy cutoffs. The errors are 
statistical and systematic. 



(GeV) 



(E) 
(GeV) 



(E 2 ) - W 
(GeV 2 ) 



1.897 
1.999 
2.094 
2.181 
2.261 



2.346±0.018 
2.338±0.010 
2.365±0.006 
2.391±0.003 
2.427±0.002 



+0.027 
0.022 
+0.020 
0.017 
+0.016 
0.010 
0.008 
0.007 
+0.006 
-0.006 



0.0211±0.0057±o:6o69 
0.0239±0.0018l°;go3o 
0.0176±0.0009io;TOie 
Q.0129±0.0003lg;ggg| 
0.0082±0.0002l°;°°g 2 



X. CONCLUSION 

We have performed a measurement of the transition 
rate of b — >• using the entire BABAR T(4S) dataset. 
We find that for E 7 > 1.9 GeV, the branching fraction is 



B(B -> X s j) = (3.29 ± 0.19 ± OA 



10" 



(8) 



We determine the pair-wise correlation between the 



where the first uncertainty is statistical and the second is 
systematic. The statistical uncertainty on this measure- 
ment is based on the sum in quadrature of the statistical 
uncertainties on each of the X s mass bin yields. This 
method of combining statistical uncertainties ensures a 
reduced spectrum dependence and is different from the 
method used in the previous BABAR sum-of-exclusives ap- 
proach where one large rrix s bin was used to determine 
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TABLE X. The correlation coefficients between (E) and the variance for the different minimum photon energies based on the 
total uncertainties (statistical and systematic). 



{E) (E 2 ) - (Ef 

(GeV) (GeV 2 ) 

E^in 1.897 1.999 2.094 2.181 2.261 1.897 1.999 2.094 2.181 2.261 



1.897 1.00 0.72 0.46 
1.999 1.00 0.71 

(E) 2.094 1.00 
(GeV) 2.181 

2.261 

L897 
1.999 

(E 2 ) - (E) 2 2.094 
(GeV 2 ) 2.181 
2.261 



0.40 0.20 -0.90 -0.66 -0.36 -0.27 -0.13 
0.65 0.35 -0.39 -0.29 -0.52 -0.40 -0.18 
0.84 0.40 -0.08 -0.25 -0.81 -0.57 -0.23 
1.00 0.67 -0.03 -0.16 -0.39 -0.48 -0.42 
1.00 0.05 0.04 0.17 0.31 -0.68 
1.00 0.51 0.13 0.10 0.00 
1.00 0.29 0.24 -0.04 
1.00 0.69 -0.10 
1.00 -0.08 
1.00 



the statistical uncertainty. This measurement supersedes 
our previous measurement using the sum-of-exclusives 
approach [To] . 

We have also measured the mean and variance of the 
photon energy spectrum. At the lowest photon energy 
cutoff (E 7 > 1.897 GeV), these values are 



(E 2 ) (Ef 



(E) = 2.346 ± O.OlSl^g GeV, 
= 0.0211 ± 0.0057^:^69 Gey2 - 



(9) 
(10) 



Finally, we have also measured the best HQET param- 
eters for two photon spectrum models. For the shape 
function models [T3] these are 



m b = 4.579 



-To^GeV/c 2 



A = 0-257^ao39 GeV 2 



(11) 
(12) 



(compared with the world averages of = 4.588 ± 
0.025 GeV/c 2 and /i 2 = 0.189^ °^ GeV 2 0), and for the 
kinetic models 1 14 these are 



m b = 4.568i°;°36 GeV/c 2 , 
/4 = 0.450 ± 0.054 GeV 2 



(13) 
(14) 



(compared with the world averages of mi, = 4.560 ± 
0.023 GeV/c 2 and & = 0.453 ± 0.036 GeV 2 0). 
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Appendix A: Correlation Coefficients for Bin- Yield 
Uncertainties 

In Table |XI| we report the correlation coefficients be- 
tween the total uncertainties reported in each mass bin. 
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